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ABSTRACT 

A  semigroup  approach  to  differential-delay  equations  is  developed  which 

seems  more  suitable  for  certain  partial  inteqro-differential  equations  than 

♦ ' 

the  standard  theory.  The  method  is  applied  to  prove  a  local-time  existence 

theorem  for  equations  of  the  form  u^  m  q(u  ^,u  ,t)  ,  where  -r-g  >  0. 

tt  xt  x  x  3u  ^ 

•  >  .  Xt 

On  a  formal  level,  it  is  demonstrated  that  the  stretchinq  of  filaments  of 

viscoelastic  liquids  can  be  described  by  an  equation  of  this  form*  t 
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SIGNIFICANCE  AND  EXPLANATION 


An  existence  theorem  for  differential  equations  of  the  form 

u*.*.  *  g(u  J_»ut,t)  with  >  0  is  proved.  Here  ufc  denotes  the 

xt  x  x  3u  .  x 

xt  t 

■history"  of  u  up  to  time  tsu^fs)  »  u^lt-s).  Equations  of  this  nature 
arise  in  the  theory  of  viscoelastic  liquids.  Specifically,  we  consider  the 


problem  of  stretching  thin  filaments.  The  constitutive  law  is  assumed  to  be 


such  that  the  stress  consists  of  the  following  three  parts:  a  pressure  term. 


a  Newtonian  contribution  and  a  "genuine"  polymeric  part,  which  is  a  smooth 


functional  of  the  strain  history.  In  a  one  dimensional  approximation  this 


leads  to  an  equation  of  the  described  form. 


The  responsibility  for  the  wording  and  views  expressed  in  this  descriptive 
summary  lies  with  MRC,  and  not  with  the  author  of  this  report. 


A  SEMIGROUP  APPROACH  TO  EQUATIONS  WITH  INFINITE  DELAY 
AND  APPLICATION  TO  A  PROBLEM  IN  VISCOELASTICITY 


M.  RENARDY 

1.  Introduction. 

Consider  a  linear  functional  differential  equation  of  the  form 
(1.D  Yt  *  Lyfc, 

where  yfc  denotes  the  history  of  y  on  a  certain  interval  [t-r,t] ;yt(s)  * 
y(t-s),  and  L  is  a  continuous  linear  operator  from  C[-r,0]  into  R. 
Equations  of  this  kind  can  be  regarded  as  evolution  problems  on  C[-r,0]  and 
can  be  studied  using  semigroup  theory.  A  generator  A  of  a  C°-semigroup  can 
be  associated  to  (1.1)  as  follows  (see  [2]): 

cN>(s ) 


(1.2) 


M>  ( s )  = 


for  -r  <  S  <  0 
for  s  =  0. 


This  approach  has  been  very  successfully  applied  to  systems  of  ordinary 
differential  equations  with  finite  delays.  If,  however,  partial  differential 
equations  and  infinite  delays  are  considered,  the  operator  A  has  some 
inconvenient  properties.  For  instance,  the  semigroup  generated  by  A  is  not 
analytic  even  for  L  =  0.  Moreover,  changes  may  occur  in  the  spectrum  of  A 
for  r  +  *.  Let  again  L  =*  0,  Then,  for  every  finite  r,  the  spectrum  of 
A  consists  only  of  the  simplt  eigenvalue  zero.  If,  however,  A  is 
considered  an  operator  on  C^m(-®,0J  (the  space  of  bounded  continuous 

D 

functions  that  have  a  limit  at  -“),  its  spectrum  is  the  whole  left  half 
plane.  Such  properties  are  disturbing,  if  one  studies  partial  integro- 
differential  equations  such  as  the  following 
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(1.3) 


ut(x,t)  «  uxx(x»t)  +  /*  a(t-s)f (u(x,s),ux(x,8))ds 

with,  e.q.,  Dirichlet  or  Neumann  boundary  conditions.  One  is  tempted  to 

reqard  (1.3)  as  a  perturbation  of  the  heat  equation  and  apply  the  results  of 

Sobolevskii  (101  to  prove  existence  theorems.  If,  however,  one  follows  the 

recipe  outlined  above,  then  even  the  heat  equation  itself  is  no  lonqer 

associated  to  an  analytic  semiqroup. 

For  these  reasons,  we  choose  an  alternative  method  of  representinq 

equations  with  infinite  delays  as  an  abstract  evolution  problem.  Iiet  F  be  a 

smooth  (possibly  nonlinear)  functional  from  C^iB(-*,0]  into  R.  Moreover, 

for  every  function  6  C^m(-“,0]  we  define  6  C^m(-°»,0]  by  T  w(s) 

b  r  b  r 

*  P  ( s  +r )  for  each  r  <  0.  Then  with  F  we  associate  the  followinq  operator 

PicJim(-»,oi  -*•  cjim(-^,0lt 


(1.4) 


F(y)(s)  -  F(T  ¥>). 

a 


There  is  some  temptation  to  rewrite  the  equation  yfc  *  F(yt)  in  the  abstract 
form  0  ■  F(y ) ,  y  6  C?'iln(-*,0].  However,  this  is  incorrect.  The  reason  is 

"  t  D 

that  functional  differential  equations  of  this  kind  are  usually  associated 
with  an  initial  condition 

(1.5)  y(t)  -  y0(t)  for  t  <  0, 

and  y0  is  not  required  to  satisfy  the  equation.  We  can,  however,  define 
f(t)  by 


(1.6) 


f(t) 


3y0(t)  t 

-  F(y^)  for  t  <  0 


TT 


for  t  >  0. 


It  will  be  assumed  that  y0  is  such  that  f  Is  continuous  (includinq  the 

point  t  ■  0),  bounded  and  such  that  lim  f(t)  exists.  Then 

t*-“ 

f(t)  6  C^m (-••,0)  can  be  defined  by 
n 


(1.7) 


f(t)(s)  -  f(t+s) 


This  permits  to  rewrite  the  equation  yt  =  P(yt )  in  the  form 
<1.8)  -  *(y>  +  t. 

In  this  paper  this  approach  will  be  applied  to  equations  of  the  form 

Utt  "  g(Uxt'Ux>x'  *  «  M'13 

<1.9) 

gtu^.u*)  -  f±1<t)  at  x  -  ±1. 
da 

It  will  be  shown  that,  if  v  a  >  0*  the  Sobolevskii  theory  of  quasilinear 

«u  . 
xt 

parabolic  equations  can  he  applied  to  this  problem  to  prove  an  existence 
theorem  locally  in  time. 

In  f3,  we  will  show  how  the  problem  of  stretching  thin  filaments  of 
polymeric  liquids  reduces  to  an  equation  of  the  type  (1.9).  For  this,  we  use 
a  formal  approximation  of  a  similar  nature  as  those  used  in  rod  theory  [7]. 
Other  problems  that  our  result  would  apply  to  are  the  exact  one-dimensional 
motions  considered  by  Coleman,  Gurtin  and  Herrera  [1]. 

Remark. 

A  little  bit  of  care  is  needed  when  proceeding  to  (1.8).  Namely,  we  must 
make  sure  that  the  abstract  solution  y(t)(s)  actually  corresponds  to  a 
function  y(t+s),  in  other  words,  Tgy(t)  ■  y(t+s).  For  this,  observe  that, 
if  y(t)  is  a  solution  of  <1.8),  then  s(t)  ■  Tgy(t)  solves  the  equation 
(1.10)  ■  fr(*l  +  Tgf. 

So  does  s’(t)  “  y(t+s).  Hence,  if  you  have  a  uniqueness  theorem  for  (1.10), 
you  are  done.  Such  a  uniqueness  theorem  will  automatically  come  out  of  the 
considerations  leading  to  our  existence  result. 
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2.  The  Main  Result 


We  study  the  following  problem 


tt 


g(u  ,ut) 
xt  x  x 


for  x  6  t— 1 ,11 


(2.1) 


g(uxt'ux5  =  f±1<t}  at  X  *  41* 

In  view  of  the  application  we  have  in  mind,  we  require  g  to  be  defined  only 
on  an  open  subset  of  the  history  space  C4im(-“,0] .  We  denote  this  open 
subset  by  Q.  We  assume  then 

(i)  g  is  a  smooth  function  from  R  x  ft  into  R.  Moreover,  there  exists 

e  >  0  such  that  D^g  >  e  on  R  x  ft.  (By  D^g,  D^g  we  denote  the 
Frechet  derivatives  of  g  w.r.  to  the  first  and  second  argument.). 
Let  us  first  study  the  boundary  conditions.  They  yield  equations  of  the 


form 

(2.2)  gtyt-.y4)  *  f(t) 

(we  have  substituted  y  for  ^(±1)).  According  to  our  assumptions  on  g, 
we  can  resolve  this  equation  w.r.  to  yt,  leading  to  an  equation  of  the 
form 

(2.3)  yt  =  h(yfc,f (t) ) 
where  h  is  smooth  from  fl  *  R  into  R. 

As  usual,  we  study  (2.2)  for  t  >  0  with  an  initial  condition  y(t)  - 
yQ(t)  for  t  <  0.  On  this  initial  condition  we  assume 

dy0(s) 

(ii)  y.  S  8  has  a  continuous  bounded  derivative,  and  lim  — t -  =  0. 

0  _  ds 

s+-°° 

Moreover,  the  compatibility  condition 
dy 

9^“dt  (0,»yo)  "  f{0) 

is  satisfied. 


-4- 


:'.v  V 


nil 


We  now  define  f(t)  by 


(2.4) 


f(t)  = 


dt  '  W 
the  given  f 


for  t  <  0 

for  t  >  0. 


As  in  the  introduction,  let  y(t)(s)  *  y(t+s),  f(t)(s)  ■  f(t+s).  Then  (2.3) 

is  rewritten  in  the  abstract  form 

(2.5)  yt  *  fi(y,f (t) ) 

where  fi  is  smooth  from  0  x  C^iin( -*•,()]  into  C^im(-",0).  The  Picard- 

D  D 

Lindelof  theorem  now  immediately  yields  the  following  result. 

Theorem  2.1. 

Let  (i),  (ii)  hold.  Then,  for  some  T  >  0,  equation  (2.5)  with  the 
initial  condition  y(0)  =  yQ  has  a  unique  solution  $  6  c\[0,T), 
c£im(--,°1). 

We  now  turn  to  the  study  of  (2.1).  It  is  convenient  to  rewrite  this 
equation  as  a  system.  We  put  u  *  p,  u  «  r,  u  =*  q.  Then  (2.1)  is 

^  t  XX 

equivalent  to 


(2.6)  q  *  r 

t  xx 

■  V'V**’  •  02«'VP*>  •  « ‘ • 

Remark. 

We  have  assumed  here  that  the  pointwise  derivative  q  =  px  is  actually 
the  derivative  in  the  sense  to  the  space  C^*"’(-<»,0).  This  is  justified  if 
q(t+s,x)  is  uniformly  bounded  and  uniformly  continuous  for  x  6  [-1,1), 
s  6  (-*•,01.  Our  assumptions  on  the  initial  history  will  he  such  that  this  is 
implied. 


-5- 


Pro*  our  study  of  the  boundary  problem,  we  know  solutions  ux(t,1)  = 
y.,(t),  ux(t,-1)  *  y_,,(t).  With  (2.6)  we  must  then  associate  the  boundary 
conditions 


(2.7) 


?x(t,D 


dy1 (t) 

~~sr~ 


?x(t,-i) 


dy-1(t) 

~k — 


dy, (t)  dy  ,(t)  ^  0,dy, (t)  dy_1(t) 

)  9 


Lettinq  r  =  r  -Vjxl — -  )  -  V4  x  ( — ^ 

we  obtain  homogeneous  Neumann  conditions  for  r . 


TT 


Before  we  reformulate  (2.6)  as  an  abstract  evolution  problem,  we  need  to 
introduce  some  function  spaces.  Let  H°  denote  the  space  of  all  real-valued 
functions  on  f-1,11,  which  have  n  square  inteqrahle  derivatives.  By 
c£im((-»,01,Hn)  we  denote  the  space  of  continuous,  bounded  functions  from 
[-•,0]  into  h",  which  converge  to  a  limit  at  We  assume 

(iii)  q  (reqarded  as  actinq  pointwise  w. r.  to  the  space  variable)  is  a 
smooth  mapping  from  R  *  0  into  H1 ,  where  ft  is  an  open  subset 
of  C3lm( (-'",01  jH1 ) .  On  R  *  ft,  we  have  D  g  >  e  >  o. 

D  I 

We  are  studying  an  initial  history  value  problem,  i.e.  u  is  qiven  for 
t  <  0.  We  assume  that  this  initial  history  is  such  that 

(iv)  ut  lies  in  C^ilB(  (--,0)  ,H3),  and  utt  lies  in  C*lm(  (--,0)  ,h’  ). 

Moreover,  u  lies  in  C^n'((-“,0),  H2  n  ft), 
x  D 

The  history  of  u  provides  histories  for  p,q  and  r  in  the  obvious  way. 

We  make  sure  that  the  history  satisfies  (2.6)  by  addinq  an  inhomoqeneous 
term  f*(t,x)  to  the  riqht  side  of  the  third  equation.  We  assume  a 
compatitility  condition 

(v)  The  initial  history  is  such  that  f*(0,x)  =  0. 

Finally,  we  assume  the  followinq  smoothness  conditions 


'  *  r-t  *-  • 


6- 


1 


f*  (regarded  as  valued  in  H1 )  is  uniformly  Holder  continuous 
in  t.  Moreover,  y,  and  y_1  (the  boundary  values  of  ux)  are 

d2yi  d2y  1 

in  C  ,  and  — ~  , - —  are  uniformly  Holder  continuous. 

dt2  dt2 


For  abbreviation,  let  us  put 


dy  (t)  dy  (t)  ,  dy  (t)  dy  (t) 

V>(t,x)  =•  V2 x(  +  dt  ^  “V4X  (  dt  dt  ^ ' 


..  .  3*(t,x)  ...  .  _  32^(t,x)  ..  .  _  »*(t,x) 

^  (t»x)  =*  9  ^2  ^,x  3  2  9  jit  * 

3x 


As  before,  we  define  £(t)(s)  »  p(t-s)  etc.  Equation  (2.6)  now  assumes  the 
following  form 

K  ■  *«  *  * , 

(2.0)  %.  .  ta  .  }2 

it  "  ^W^xx^1  +  °2^(?x+^1  ”^3  +  f 


with  the  boundary  condition  »  0  at  x  *  ±1. 

Me  quote  the  following  result  of  Sobolevskii  (10) 
Theorem  2.2. 


Consider  an  evolution  equation 


(2.9) 


vt  +  A(  trv)v  =  f(t,v) 


in  a  Banach  space  E  .  Assume  the  following  are  trues 

Aq  «  A(0,Vg)  is  a  linear  operator  whose  domain  D  is  dense  in  E. 

Let  there  be  d  6  *  such  that,  for  Re  A  >  d,  the  operator  AQ  +  Xl  has  a 
bounded  inverse,  and 


(2.10) 


ma0  +  xi)-1 1  <  c( |x|  + 1 r1. 

For  some  a  6  (0,1)  and  for  anv  v  6  E,  Ivl  <  R,  let  the  operator 
-a 

A(t,AQ  v)  he  defined  in  D  and 

(2.11)  I [A( t,A^°v )  -  A(T ,a”°w) Ja”1 I  <  C(R)[|t-r|e  +  Iv-wl),  e  6  (0,1) 
for  anv  0  <  t,T  <  T,  Ivl  <  R,  Iwl  <  R.  in  this  reqion,  let 

(2.12)  If (t,A~°v)  -  f (T,A~°w) »  <  C( R ) ( | t-T | 6  +  Iv-wl). 

B  ci 

Finally,  for  some  8  >  o,  let  vQ  6  D(Ag),  and  let  lA^Vgl  <  R*  Then  there 
exists  one  and  only  one  solution  of  equation  (2.9)  which  is  defined  on  an 
interval  (0,tQ),  is  continuous  for  t  6  [Ojtg),  is  continuously  differ¬ 
entiable  for  t  >  0,  and  satisfies  the  initial  condition 

(2.13)  v(0)  =  vQ. 

We  apply  this  to  (2.8)  by  makinq  the  following  identifications.  Let  E 
denote  the  space  C*im( (-**,0) ,H* )3  and  let  v  =  (p-p  ,q-q  ,r-r  ),  where 

n  0  0  0 

P0,<50,*o  denote  the  qiven  initial  conditions.  The  operator  A(t,v)  is 
defined  by  -A( t, v) (p* ,q ' ,r ' )  =  (?',?•  ,D  q(r  +£.,$)£’  )  incorporating  the 

X  XX  I  X  I  XX 

boundary  condition  r^  =  0  and  f(t,v)  is  identified  with  the  remainder  of 

the  right  hand  side  of  (2.8).  Then  the  initial  condition  vQ  =  0  is  in 

D(A0),  i.e.  we  can  choose  8=1.  Moreover,  for  appropriate  choice  of  R,  it 

is  clear  that  (2.11),  (2.12)  hold  with  o  =  V2  •  and  e  the  Holder  exponent 

of  condition  (vi).  Condition  (2.10)  finally  follows  from  the  positivity  of 

D^q  and  the  fact  that,  for  positive  ^(x),  the  operator  r  ►*  ^(x)rxx  with 

homogeneous  Newmann  conditions  is  selfad joint  in  H1,  provided  that  H1  is 

equipped  with  the  scalar  product  ((p,q))  = 

J1  — J — r  p(x)q(x)dx  +  /'p  (x)q  (x)dx.  This  product  is  equivalent  to  the  usual 
.  t  1 X I  ,  X  X 


one 


We  thus  obtain 


Theorem  2.3. 

Let  hypotheses  (i)-(vi)  hold.  Then,  for  some  T  >  0,  equation  (2.8)  has 
a  unique  solution  in  C^IO.T],  C^*m(  -®,  01  ,H^ )  3 ) ,  which  satisfies  the  qiven 

D 

initial  conditions. 

Remarks . 

Certain  models  for  viscoelastic  liquids  express  the  stress  as  a 
functional  not  only  of  the  history  of  strain,  but  also  of  the  history  of  the 
rate  of  strain  (see  e.q.  111]).  In  this  case  we  will  have  to  deal  with 
equations  of  the  form 

(2.4)  utt  =  <3<uxt' VUxt)X* 

If  D^g  is  small  compared  to  D.q,  it  can  be  shown  by  perturbation  arguments 
that  the  assumptions  necessary  for  the  theorems  above  remain  still  valid. 

In  many  cases  one  has  to  deal  with  functionals  that,  althouqh  they  involve  the 
rate  of  strain,  also  have  a  smoothinq  property.  The  stress  is  (at  least)  as 
often  differentiable  w. r.  to  time  as  the  strain,  and  no  derivative  is  lost. 
E.g.  take  a  constitutive  law  of  the  form  a(t)  *  /^(t-sjffu  (s))ds.  This 

_oo  xt 

leads  to  a  (t)  =  /*  a'(t-s)f(u  i_(s))ds  +  a(0)f(u  ft)).  Thus  o  is  given 

t  '  xt  xt  t 

by  an  expression  involvinq  only  uxfc  and  not  uxtt*  the  h*story  uxt 
enters  (2.4)  via  such  expressions,  it  is  possible  to  apply  the  same  methods  as 
before  to  the  time  differentiated  version  of  (2.4).  The  terms  involving 
histories  are  then  still  of  "lower  order"  and  can  be  included  in  the  "f- 
part".  We  leave  the  details  to  the  reader. 
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3.  The  Stretching  of  Thin  Filaments  of  Viscoelastic  Liquids. 

We  consider  the  problem  of  stretching  a  thin  filament  of  a  viscoelastic 
liquid  by  a  force  f  applied  to  its  ends: 


We  shall  derive  an  equation  describing  this  problem  in  an  analogous  fashion 

as  I  did  in  18]  assuming  a  special  constitutive  law.  We  start  with  the 

general  three-dimensional  equations  of  continuum  mechanics  and  then  introduce 

a  one-dimensional  approximation  analogous  to  rod  theory  [7). 

12  3 

Let  (C  ,C  ,5  )  denote  "body  coordinates",  i.e.  coordinates  labelling  a 

specific  particle  in  the  fluid.  These  coordinates  can  be  identified  with  the 

position  of  the  particle  in  space,  when  the  fluid  is  in  a  certain  "reference 

state".  This  reference  state  will  be  chosen  such  that  the  filament  is  cylindrical 

in  this  state.  By  (y  ,y  ,y  )  we  denote  Cartesian  coordinates  in  space.  We 

We  are  interested  in  finding  trajectories  of  particles,  i.e.  a  functional 
12  3 

dependence  y :  ( £  ,G  ,t). 


The  deformation  gradient  F 


i  dy 

is  defined  by  F  j  = 

3C3 


and  the  right 


Cauchy-Green  strain  tensor  is  defined  by 


(3.1) 


C  =  FTF. 


Following  Lodge  [3),  we  introduce  a  Lagrangian  stress  tensor  S  ,  which  is 
related  to  the  Cauchy  stress  tensor  T  by 

(3.2)  s  -  F-Wr1. 

S  is  related  to  C  by  a  constitutive  law.  We  assume  that  this  constitutive 
equation  consists  of  three  terms:  a  pressure  term,  a  Newtonian  viscosity  and 
a  "genuine"  polymeric  contribution,  which  is  a  smooth  functional  of  the  strain 
history.  T.e.,  we  have 

(3.3)  S  .  -p  c'1  -  n  c_1  +  F(Ct). 
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The  incompressibility  condition  is  expressed  by 

(3.4)  det  C  *  1. 

The  equation  of  motion  is  then  given  by 

(3.5)  pytt  -  Div(FS) 

where  Div  refers  to  the  body  coordinate  system,  or,  in  coordinates 


(3.6) 


-I 

r,s 


.  ,2  It 

_J_  s«  +  _1jl_  Srs. 

3CS  3 Cr3?8 


We  have  to  supplement  (3.6)  by  the  boundary  conditions.  Let  f  denote  the 
force  acting  on  a  surface  divided  by  the  surface  area  in  the  reference 
state,  f  is  referred  to  space  coordinates.  Then  the  boundary  condition 
reads 

(3.7)  P  •  S  •  n  -  f, 

where  n  denotes  the  unit  normal  of  the  surface  in  the  body  coordinate 
system. 

We  now  specialize  our  considerations  to  the  filament.  Coordinates  are 
chosen  such  that  C1  and  y1  are  in  the  direction  of  the  filament,  and  the 


other  coordinates  are  transverse*  It  is  assumed  that  in  the  reference  state 

the  filament  is  axisymmetric  and  cylindrical,  i.e.  51  ranges  over  a  finite 

/  2  2  3  2 

(C  )  +  (C  )  ranges  from  0  to  6,  the 

radius  of  the  filament.  We  shall  assume  that  6  is  small.  In  a  formal 
argument,  we  shall  derive  an  equation  for  the  terms  of  highest  order  w. r.  to 
6.  Hie  argument  is  analogous  to  the  one  given  by  Nariboli  J7)  for  the 
problem  of  longitudinal  elastic  waves  in  a  thin  rod. 

We  assume  that  the  solutions  y1,y2,y3  and  the  pressure  p  are  such 


that  they  depend  smoothly  on  r,  and  that  the  coefficients  in  a  Taylor 


expansion  w.r. 


to  powers  of  r  are  smooth  functions  of  6.  Let  us  put 
3 

■j— ,  ?  ■  so  that  the  lateral  surface  corresponds  to 
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?=1 .  For  the  solutions  we  now  make  the  Ansatz: 

1 .;3,  r  ►  2V  1  a2.  v,r2N+2. 

y  (C  ,6?  ,o?  )  =  l  5  P  (C  ,r  )  +  0(6  ) 

v=0 


N 


=5^  l  62\{C\r2)  +  <*«2N+3) 

v=0 


(3.8) 


3 .  1  .  2  j_3.  f-3  r  . 2V  1  A 2 .  ^.,2N+3. 

y  (C  ,«C  /<5?  )  =  2.  5  Qv(c  #r  )  +  0(6  ) 

\>=0 

P<c\«C2,S?3)  =  l  62VR  (c1,^2)  +  o'(«2N+2) 

V=0 

where  Pv,(?v,Rv  are  polynomials  of  vth  degree  in  r2.  We  insert  this 
Ansatz  into  (3.4),  (3.6)  and  the  boundary  conditions  (3.7)  on  the  lateral 
surface  of  the  filament.  By  retaining  only  lowest  order  terms  in  6,  we 
shall  obtain  an  equation  for  Pg(C^).  In  the  following  analysis  we  shall 
carry  out  the  series  expansions  (3.8)  only  so  far  as  is  necessary  for  this 
purpose. 

For  the  Cauchy-Green  tensor  we  find 


Thus,  in  the  lowest  order,  the  incompressibility  condition  (3.4)  yields 


(3.10) 


,*V2  4 

3  5  0 


1. 


Next,  let  us  consider  the  boundary  conditions  on  the  lateral  surface.  At  a 


boundary  point  where  5  =0  these  yield  the  equations  S21  =  S22 


.23 


•  ***.*<' 
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(Because  of  the  radial  symmetry  it  suffices  to  consider  these  boundary  points; 


if  the  traction  vanishes  there,  it  does  so  everywhere  else),  S2^  vanishes 
identically  as  a  result  of  the  radial  symmetry.  For  S  we  obtain  the 
following  terms  of  order  0(1): 

(3.11)  s22  =  -R0Q-02  -  n^(Q-2)  +  F22(cJ) 

where  CQ  denotes  the  zeroth  order  term  in  C.  Moreover,  we  have  to  satisfy 

(3.12)  S21  *  0. 

It  follows  from  the  principle  of  material  frame  indifference  that,  if  C  is 
diagonal,  then  S  must  he  diagonal,  too  (see  e.g.  [3)).  Hence  (3.12)  is 
satisfied  in  0(6),  if  f  =  0.  The  equation  of  motion  now  yields  the 
following  equation  for  P^: 


(3.13) 


3p„  »  3P^ 

0/3  ,  _  f  0  i-2 

i  i  -v  i 

3C1  3c  3? 


32p„ 


ocV 


(co),} 


+  F  (c0))« 


In  order  to  simplify  notation,  let  us  write  u  for  PQ  and  x  for  C  . 

Since  Cn  depends  only  on  u  ,  we  shall  also  write  F1 1  (c!!)  ^F^tuS, 

22  t  2  t 

F  (CQ)  =  F  (u^).  When  we  substitute  Rg  from  (3.11),  we  obtain  the 

following  equation  from  (3.13) 

'2  ,  1  ,  3  ,  c1,  t.  -2C2,  t.x 

( - )  +  ■*— iu  F  (u  )  -  u  F  (u  )}. 


(3.14) 


PUtt  =  ^SxSt 


fx  X 


X 

Finally,  we  have  to  derive  the  boundary  conditions  at  the  ends  of  the 
filament.  As  noted  in  [7],  expansions  of  the  type  (3.8)  can  not  be  used  to 
satisfy  these  in  every  ordec  in  6,  and  a  "boundary  layer"  arises.  However, 
for  the  first  order  approximation,  we  need  not  ileal  with  this.  When  we  insert 
(3.8)  into  the  boundary  conditions  at  the  ends,  then  all  traction  components 
not  in  the  direction  of  the  Filament  are  <*(6).  For  the  longitudinal 


traction ,  we  find  in  order  0(1): 

(3.15)  3H  tt-  (-  — )  +  tu  F1{ut)  -  u“2F2(ut)}  «=  f(t) 

«t  u  x  x  x  x 

x 

where  f(t)  denotes  the  force  acting  on  the  end,  divided  by  the  cross- 

2 

sectional  area  in  the  reference  state,  *5  . 

Remark. 

Our  analysis  in  chapter  2  breaks  down  if  n  *  0.  However,  we  have  seen 
from  numerical  studies  15 J,  (61  of  a  specific  model  that,  for  certain  fluids 
(like  molten  polyethylene),  >1  must  be  assumed  very  larqe  in  order  to  have 
significant  influence  on  the  solution.  Experimental  values  for  n  are  not 
available.  Thus  a  theory  that  includes  the  case  h  =*  0  seems  to  have  more 
than  academic  interest.  The  equations  are  not  parabolic  in  this  case,  but  can 
for  a  wide  class  of  constitutive  laws  be  classified  as  hyperbolic.  This  will 
be  discussed  in  a  forthcoming  paper  (9).  In  a  restricted  situation.  Lodge, 
McLeod  and  Nohel  (41  have  discussed  solutions  to  the  boundary  equation  (3.15) 
in  the  limit  n  ♦  0. 
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